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1 We compute the difference in decay widths of charged and neutral p(770) vector 

mesons. The isospin breaking arising from mass differences of neutral and charged 
7r and p mesons, radiative corrections to p — > tttt, and the p — > ttitj decays are taken 
Q H i into account. It is found that the width difference Ar p is very sensitive ot the isospin 

breaking in the p meson mass Am p . This result can be useful to test the correlations 
observed between the values of these parameters extracted from experimental data. 
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Lowest-lying vector mesons undergo predominantly strong interaction decays. The masses 
and decay widths of members of the same isomultiplet will therefore look very similar [1| , with 
small differences induced by the breaking of isospin symmetry. The isospin breaking effects 
in the p meson parameters have raised an interest recently, due to both experimental and 
theoretical reasons . According to the PDG 

measurements are: 



X 



1] the weighted averages of available 



Am p = m p o - m p ± = (-0.7 ± 0.8) MeV , (1) 
AT p = F> - T p ± = (0.3 ± 1.3) MeV . (2) 

These results are consistent with the absence of isospin breaking in the p° — system. Note 
lowever that the scale factors associated with the above averages are, respectively, 1.5 and 1.4 
l| which reflects an important spread in the yields from different experiments. 

Some recent theoretical calculations of Am p seem to confirm the above result. Using a 
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vector-meson dominance model to parametrize the jpp vertex, the authors of Ref. [4j have 
obtained Am p = (—0.02 ± 0.02) MeV. Also, using 1/N C expansion techniques, the authors 
of Ref. [3|] have obtained —0.4 MeV < Am p < 0.7 MeV. On another hand, it has been 
found that the width difference of p mesons is of great importance to understand the current 
discrepancy between the hadronic vacuum polarization contributions to the muon anomalous 
magnetic moment obtained from r decay and e + e~ annihilation data j^, 5, 6]. 

In this paper we provide an estimate of Ar p by considering the isospin breaking corrections 
in the exclusive modes that contribute to the decay widths of p 0,± vector mesons. A previous 
estimate of this effect was done in reference taking into account several sources of isospin 
breaking as mass differences and other subleading p meson decays. Their result Ar p « (— 0.42± 
0.59) MeV I2I is consistent with the world average given in Eq. (2). Additional contributions 
to isospin breaking in Ar p , including the radiative corrections to the dominant p — > tttt decays, 
are considered in this paper. 

II. SOURCES OF ISOSPIN BREAKING 

At a fundamental level, isospin symmetry is broken by the different masses of u and d quarks 
and by the effects of electromagnetic interactions. At the hadronic level all manifestations of 
isospin breaking can be traced back to such fundamental sources. In the absence of isospin 
breaking, the p 0,± mesons must have equal masses and decay widths, thus Am p = AT p = 0. 

The dominant decay modes of p mesons that are common to charged and neutral p's are the 
mi decay and its radiative mode. The branching fraction of other modes contributing only to 
the p° meson adds up to 

B Q rest = 5(7T% 7 ) + £(777) + B(fi + p-) + B(e + e-) + 5(7r + 7r-7r ) w 5.3 x 10" 4 . (3) 

There is also a dipole transition p — * tij which is common to p ±,Q vector mesons with branching 



fractions of a few times 10 4 



Since the p meson widths are of order 150 MeV, all these 



subleading decay modes will contribute to the width difference at the tiny level of: 

Ar sub w 0.O8 MeV . (4) 
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Thus, any sizable difference in the decay widths can only originate from the dominant decay 
modes. To be more precise, we will define explicitly the contributions to the width difference 
as follows: 

Ar p = r(p° 7r + 7r-( 7 ),^ < LU ) - T(p + -> 7r + 7r°( 7 ),^ < oo ) 

+ r(p° -> 7T + 7T-7, UO > U0 ) ~ T{ P + -> 7T + A, ^ > ^o) + AI^ . (5) 

The first two terms in Eq. (jSJ) denote the tttt decay rates that include virtual and soft-photon 
corrections and the next two terms are related to hard-photons in TTTc-f decays. In the above 
equation, uu is an arbitrary value of the photon energy that separates the decay rates of soft- 
and hard-photon bremsstrahlung. It is expected that the ujq dependence will cancel in the sum 
of the first (second) and third (fourth) terms of Eq. (jSJ). 

In the following section we will consider in more detail the isospin breaking corrections to 
p — > 7T7T and their radiative decays. Later, we will evaluate the contributions of such corrections 
to the r.h.s. of Eq. (JSJ). 



III. RADIATIVE CORRECTIONS TO p -> vrvr DECAYS 
At the lowest order (indicated with superscript 0), the rates of p — ► titt decays are given by: 

r° = rV - TT+O = f±-mX , (6) 

I* = r°(p+ - vr+vr ) = ^m p+ vl , (7) 

4o7T 



where 



4m 2 + 

Vq = \ 1 1 1— , and f + = 



\ 



^ (m 7r + — m^o) 2 ^ ^ (m T + + m^o) 2 



2 II 2 

m i / \ m 



are the pion velocities in the rest frame of the p l meson and are the p7rV J coupling constants 
such that <7_|__ = (7+0 owing to the isospin symmetry of strong interactions. In the limit of 
isospin symmetry, the masses of neutral and charged pions (p's) are the same, v = v + , and 
consequently the three-level decay rates are equal (rfj = TV). In the following subsections we 
will consider the different ingredients to get the 0(a) radiative corrections to the decay rates 
given above. 
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FIG. 1: Feynman graphs for p° — > 7r + 7r 7 decays. 



A. Real soft-photon corrections 



In order to get an infrared safe result, the radiative corrections of order a to Eqs. (6,7) must 
include the sum of virtual corrections and soft-photon bremsstrahlung. In order to define the 
soft-photon contribution, let us consider the radiative p — > tttc-/ decays shown in Figures 1 and 
2. In each of Figures 1 and 2, the first three diagrams correspond to the model- independent 
contributions and the other graphs denote the model-dependent terms. 

The decay amplitude of the process p°(d,rj) — > it + {p)it~ {p'^ik, e*) (four-momenta and 
polarization four- vectors are indicated within parenthesis), see Figure 1, is given by [7, 8]: 

p ■ e*k- 1] 



(p° — > 7r + 7r 7) = ieg + - < (p — p ) ■ 7] 



p ■ e 



p ■ e 



+ 2 



e ■ 77 — 



p ■ k 



+ M d , e>i (9) 



where the model-independent pieces of the amplitude are given explicitly, and the remaining 
terms -Md,e,f correspond to the model-dependent amplitudes arising from Figures le)-lf). Sim- 
ilarly, the decay amplitude for the decay p + (d, rj) — > 7r + (j9)7r°(p / )7(A;, e*)), see Figure 2, is given 
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FIG. 2: Feynman graphs for p + — > 7r + 7r°7 decays 
by the following expression [7, 9, 1C 
.M(p + — ► 7r + 7r°7) = ie^+o 



111: 



p.e 



d.e* \ ( p.e 

{p — p)-r) + 



+ 



p.k d.k 



2 + ^d + ^ 



d.e* 



m 



d.A; 



■k ■ 7] — e* • r] 



k ■ i] 



-(2 + An) [ —k ■ r] — e ■ rj 



p ■ k 



d ■ k 



M 



d,e ) 



(10) 



where Ak denotes the anomalous magnetic dipole moment of the p + vector meson in units of 
e/2m p + (in our numerical evaluations we will set Ak = 0), and we have defined = 



-m 2 . 



Once again, only the model-independent amplitudes for Figures 2a)-2c) have been written 
explicitly. 

The soft-photon bremsstrahlung is the divergent piece of the amplitudes (jnHHD in the infrared 



region {k — > 0). The diagrams contributing to the soft-photon amplitude are shown in Figs. 
la)-lb) for p° — > 7r + 7r~7 and in Figs. 2a)-2b) for p + — > 7r + 7T°7 decays. As usual, only photons 
of energy smaller than uq must be considered in the real-photon emission rate that is included 
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in radiative corrections. Thus, the radiative corrected rates of p l — ► tttt decays (i = 0, +) can 
be written in the following general form: 



T(p l -> 7T7r(7),u; < lUq) 



r o + r , lv + rr ft M 
r? 



r° 







where r| v denotes the virtual corrected rate at order a for p* — > 7T7t decay and (a; ) is the 
soft-photon rate of p % — > 7T7T7 obtained from the first term in Eqs. (MTU]) . 

The soft-photon rates can be computed in analytical form by integrating over photon energies 
smaller than ujq, with the following results, respectively for p° and p + radiative decays: 



"psoft 
1 



r° 

1 



7T [ \2lU 



1 + 



1+1* 



2^n 



hi 







2Vn 



Li 2 



^0 



+ ln 



l + ^o 
In 



Li 2 



l + v 
1+vp 

1-Vq 



iln 

^0 



+ Z7rln 



1 + ^0 



(12) 



and 



"'soft . 



r° 




(13) 



where we have defined v' + = i>+/6> with £ = l + A^/m^ + , and Li(x) denotes the Spence function. 

Observe that the soft-photon rates depend logarithmically upon the infrared cutoff A (the 
photon mass regulator) and the photon energy parameter ujq. We have checked that Eq. (12) 
coincides numerically with the result reported long ago by Cremmer and Gourdin [12| in the 
case of <p —>■ K + K~ r y decays, although our corresponding expressions are written in a different 
form. Finally, let us point out that a very small contribution arising from the regular part 
(finite when the photon energy goes to zero) of the radiative decay rate must be added to the 
r.h.s. of Eq. (11). 
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FIG. 3: Feynman graphs for the virtual photonic corrections to p° — > ir + ir decays 



B. Virtual corrections 



The Feynman diagrams of virtual corrections of 0(a) to the p —>■ nix amplitudes are shown 
in Figures 3 and 4. As in the real photon case, we will use scalar QED for the electromagnetic 
couplings of charged pions and assume that the p electromagnetic vertex (see for example 
0,0]) is similar to the WW'y gauge boson vertex. 

The sum of virtual corrections to the p° — > ix + ix~ amplitude is finite in the ultraviolet region 
owing to a Ward identity satisfied by the vertex and self-energy corrections 121]. In the case 
of the p + — > 7T + 7r° decay, the sum of virtual corrections diverges. One way of getting rid of 
such divergences is by absorbing them into a redefinition of the strong coupling constant g + Q. 
Instead of dealing with divergent amplitudes, in the case of p + — > 7r + 7r° decays we will follow 
a procedure introduced long ago by Yennie 13j. According to this method we pick only the 
convection terms from virtual corrections in Figures 4a)-c). The sum of such terms gives rise 
to a one-loop amplitude which is finite in the ultraviolet region and is also gauge-invariant and 
model-independent jl3]. Furthermore, such an amplitude contains all the infrared divergent 
terms of virtual corrections [3]. 
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FIG. 4: Feynman graphs for the virtual photonic corrections to p + — > tt + tt° decays. 



Just to illustrate the cancellation of infrared divergencies in the sum of virtual and real 
photon corrections, we reproduce here the expression for the virtual corrections in the case of 



P 



7T + 7r decays 12]: 



"plv 
1 

r° 

1 



a 

7T 



71 



2Vn 



2 1 + In 



A 



77V+ 



1 + 



2Vn 



vo 



[Li 2 (v ) - Li 2 (-t) )] - 



2Vn 





2 


(u 2 


l + v _ 



l+v 
- Li 2 



• (14) 



The corresponding result in the case of p H 
described above [13] . We obtain: 



7r + 7T° decays can be computed using the procedure 




7T 

— + 2Li 2 . 
3 V 1 + v' 
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^(^) + -( zl± Sr M )]} 



(15) 



where we have defined = m^ + +m^ - As pointed out before, this result is free from ultraviolet 
divergencies, and it depends logarithmically on the infrared mass regulator A. 

As we can easily check, the sum of the virtual corrections, Eqs. ( ll4|[T5"j) and the corresponding 
soft-photon bremsstrahlung, Eqs. (I12lll3p . are free from the infrared regulator A as it must be 
according to the Bloch-Nordsiek theorem. Note however, that this sum depends on the photon 
energy cut uq- We can expect that this ujq dependence will be largely cancelled in the photon- 
inclusive p — ► 7T7T decay rate in such a way that Ar p defined in Eq. (5) does not depended on 
this parameter. Since the soft- and hard-photon decay rates have a different dependence upon 
ujq, we need to choose values of ojn which are sufficiently small (typically much smaller than the 
masses of hadronic particles) [14|. In the following subsection we briefly comment about the 
decay rates for hard photons. 



C. Hard-photon emission 

In order to estimate the contributions to Ar p due to photons of energy larger that ujq, the 
third and fourth terms in Eq. (5), we require the complete expressions of the radiative decay 
amplitudes corresponding to the Feynman diagrams in Figures 1 and 2. These amplitudes 
include also the model-dependent contributions shown in Figures ld)-lf) and 2d)-2e). The 
explicit contributions of the model-dependent amplitudes can be found in Refs. [3, Rill |- 

When we evaluate the decay rates of radiative decays, we note that the contributions of 
the model-dependent terms (intermediate states with u>, ai,a and f mesons) are typi cally 
two orders of magnitude below the contributions due to model-independent terms jsl, O, \lu\ . 
Therefore, such model-dependent terms would affect the width difference Ar p at a completely 
negligible level. Note however that their effects have been included in our numerical evaluations 
given below. 
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IV. RESULTS AND DISCUSSION 

In this section we present the main numerical results of our calculations. The values of 
the radiative corrections 5 p i (i = +, 0) defined in Eq. (11) are shown in Tables I and II for 
p° — > 7r + 7r~(7) and p + — > 7r + 7r°(7) decays, respectively. The values of 5 p i are tabulated as 
a function of the photon energy cut c^o for three different values of the p % meson mass which 
are consistent with a small isospin breaking. As expected, these radiative corrections depends 
only slightly on the specific value of the p meson mass. It is interesting to observe that in 
cf) — > K + K~ decays (and other analogous quarkonia decays that occur close to threshold) the 
radiative corrections are dominated by the coulombic interaction in the final state [12|], while 
in p° — > 7r + 7r~ this is not the case. 

Table III shows the radiative decay rates of p l — > tctc^ normalized to the the corresponding 
tree-level rates The tabulated values are defined as follows (i = 0, +): 

_ IV -» 7T7T7,^ > Up) 

A p< = fo • ( 16 ) 

i 

These values were calculated for photon energies larger that uq and for three different values 
of the p meson mass. The values of A p i are always positive and exhibit the typical decreasing 
behavior as the lower energy photon cut ujq increases. We observe that these normalized rates 
are in good agreement with the corresponding results reported long ago in Ref. using slightly 
different values for the mass and width of the p meson. 

The predicted branching fraction for the neutral radiative mode: 

B(p° -> 7r + 7r"7,w > 50 MeV) = 11.5 x KT 3 , (17) 

compares very well (they are in agreement within la) with the experimental measurement 
(9.9±1.6) x 10~ 3 reported in [II], which was obtained for the same value of co>o- For the same value 
of u , we obtain the isospin breaking in the radiative modes to be (assuming T p+ ps T p o = 150 
MeV): 

T(p° -> 7r + 7r"7, lo = 50 MeV) - T(p+ -> 7r + 7r°7, tu = 50 MeV) w 1.1 MeV , (18) 

which is larger that the estimate (0.45 ± 0.45) MeV considered in Ref. [s| and essentially 
independent of the p meson mass. 
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In Table IV we display the sum of the radiative corrections 5 pl and the radiative branching 
ratios A p i for a common value (775 MeV) of the neutral and charged p meson mass. According 
to the definitions given in the previous sections we have: 

T(p l -> 7T7r(7), UJ < LUp) + T(p l -> 7T7T7, U) > LUp) 

-g — i + o p i + lApi 

i 

= 1 + a p i . (19) 

As we have pointed out in section II, we expect that the ay correction term will be independent 
of the photon energy cut Up as far as Eq. ( fl~9l) describes a photon inclusive rate and up is an 
arbitrary reference value used to separate soft- and hard-photon emission. In practice this 
cancellation is better for values of Up (typically below a few MeV's, see Table IV) because 
radiative corrections includes only the logarithmic dependence in up, while the hard-photon 
radiative rate contains also linear and other ujp dependent terms. 

Finally, we can evaluate the difference in decay widths arising from different corrections. 
Based on previous definitions, we can write the width difference of rho mesons, Eq. (5), as 
follows: 

r / 3 \ ~\ 

f AF ub . (20) 



ay — r° 



(Tfl + \ 
m p0 Vp> J 



Once the isospin breaking in Am p is known, the values of Ar p can be easily evaluated from 
Tables I- 1 1 1 , using the above expression. For illustrative purposes, we provide the values of Ar p 
for two interesting cases (we have used Tq = 150 MeV and the value of a p i at Up = 10 MeV): 



Ar p 



0.86 MeV, if Am p = 



0.02 MeV, if Am p = -3 MeV . (21) 
1.70 MeV, if Am p = +3 MeV 



Note that the width difference is very sensitive to the size and sign of Am p . This is very 
interesting because it can be useful as a test of (Am p (exp)), Ar p (exp)) correlated values 
extracted from fits to experimental data. In particular, the above results are in good agreement 
with the values (—2.4 ± 0.8, —0.2 ± 1.0) MeV extracted from a combined fit to r and e + e~ 
data [la ], but the agreement is less good with the central values (—3.1 ± 0.9, —2.3 ± 1.6) MeV 
reported in {(J. 
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V. CONCLUSIONS 

The isospin breaking in the p° — p^ system is a very important ingredient to understand the 
current discrepancy in predictions of the hadronic vacuum poplarization contributions to the 
muon magnetic moment based in r decay and e + e~ annihilation data. In this paper we have 
evaluated the difference in decay widths (Ar p ) of p(770) vector mesons. We have considered 
in our calculation the isospin breaking in the exclusive decay modes of charged and neutral 
p mesons. In particular, we have carried out a calculation of the radiative corrections to the 
dominant p —>■ decay modes and have done a careful reevaluation of the differences in their 
radiative p — > mcy decays. 

We found that Ar p is sensitive to the isospin breaking in the p° — p^ 1 mass difference. This 
provides a useful tool to test the (model-dependent) values of these parameters extracted from 
experimental data. In particular we have found that positive values of Ar p are favored by 
|Am p | < 3 MeV. 
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m p0 = 772 MeV 


m p o = 775 MeV 


m p o = 778 MeV 


uj (MeV) 


S p o 


S p o 


5 p o 


2 


-0.03670 


-0.03692 


-0.03714 


4 


-0.02910 


-0.02930 


-0.02949 


6 


-0.02465 


-0.02483 


-0.02501 


8 


-0.02150 


-0.02167 


-0.02183 


10 


-0.01905 


-0.01921 


-0.01937 


12 


-0.01705 


-0.01720 


-0.01736 


14 


-0.01536 


-0.01550 


-0.01565 


16 


-0.01389 


-0.01403 


-0.01477 


18 


-0.01260 


-0.01273 


-0.01287 


20 


-0.01144 


-0.01157 


-0.01170 


30 


-0.00697 


-0.00708 


-0.00720 


40 


-0.00378 


-0.00388 


-0.00399 


50 


-0.00130 


-0.00139 


-0.00150 


60 


0.00074 


-0.00065 


0.00056 


70 


0.00249 


0.00240 


0.00232 


80 


0.00401 


0.00393 


0.00384 


90 


0.00536 


0.00529 


0.00521 


100 


0.00659 


0.00651 


0.00643 



TABLE I: Radiative correction <Lo to the p° — > it + it (7) decay rate (see definition in Eq. ([5])) as a 
function of loq and for three different values of m p o . 
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m p+ = 772 MeV m p - 


h = 775 MeV 


m p + = 778 MeV 


cj (MeV) 


6 p + 


6 p + 


5 p + 


2 


-0.01959 


-0.01968 


-0.01970 


4 


-0.01701 


-0.01710 


-0.01718 


6 


-0.01551 


-0.01558 


-0.01566 


8 


-0.01444 


-0.01451 


-0.01459 


10 


-0.01361 


-0.01368 


-0.01375 


12 


-0.01293 


-0.01300 


-0.01307 


14 


-0.01236 


-0.01242 


-0.01249 


16 


-0.01186 


-0.01192 


-0.01199 


18 


-0.01142 


-0.01149 


-0.01155 


20 


-0.01103 


-0.01109 


-0.01115 


30 


-0.00953 


-0.00958 


-0.00963 


40 


-0.00844 


-0.00849 


-0.00854 


50 


-0.00761 


-0.00765 


-0.00769 


60 


-0.00692 


-0.00696 


-0.00700 


70 


-0.00633 


-0.00637 


-0.00639 


80 


-0.00582 


-0.00584 


-0.00589 


90 


-0.00536 


-0.00540 


-0.00544 


100 


-0.00495 


-0.00499 


-0.00502 



TABLE II: Radiative corrections 5 p + to the p + — > 7r + 7r°(7) decay rate (see definition in Eq. ([5])) as 
a function of u>q and for three different values of m p % . 
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m p +,o = 772 MeV 




m p +,o = 775 MeV 




m p+t0 = 778 MeV 




oj (MeV) 


V 




Ap + 


Apo 


Ap + 


Apo 


2 


0.01544 


0.04475 


0.01553 


0.04497 


0.01561 


0.04518 


4 


0.01290 


0.03724 


0.01297 


0.03742 


0.01305 


0.03761 


6 


0.01143 


0.03288 


0.01149 


0.03305 


0.011556 


0.03322 


8 


0.01039 


0.02981 


0.01045 


0.02997 


0.01051 


0.03013 


10 


0.00959 


0.02745 


0.00965 


0.02760 


0.00970 


0.02775 


12 


0.00894 


0.02553 


0.00900 


0.02568 


0.00905 


0.02582 


14 


0.00840 


0.02393 


0.00845 


0.02406 


0.00850 


0.02420 


16 


0.00793 


0.02255 


0.00798 


0.02268 


0.00803 


0.02281 


18 


0.00753 


0.02134 


0.00758 


0.02147 


0.00762 


0.02159 


20 


0.00717 


0.02027 


0.00721 


0.02039 


0.00726 


0.02051 


30 


0.00581 


0.01624 


0.00585 


0.01635 


0.00589 


0.01645 


40 


0.00488 


0.01350 


0.00492 


0.01359 


0.00495 


0.01369 


50 


0.00420 


0.01146 


0.00423 


0.01155 


0.00426 


0.01163 


60 


0.00366 


0.00987 


0.00369 


0.00994 


0.00372 


0.01002 


70 


0.00322 


0.00857 


0.00325 


0.00864 


0.00327 


0.00871 


80 


0.00286 


0.00750 


0.00288 


0.00757 


0.00291 


0.00763 


90 


0.00255 


0.00659 


0.00257 


0.00665 


0.00259 


0.00672 


100 


0.00228 


0.00582 


0.00230 


0.00588 


0.00232 


0.00593 



TABLE III: Decay rates of p % — > ttttj (normalized to the tree-level rates F®, see definition in Eq. (16)) 
as a function of the low photon energy cut u>q and for three different values of m p . 
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lo (MeV) 








Op+ 


(TpO - a p + 


2 


8.05 x 


10" 


-3 


—A 1 ^ * 


1 0" 

1U 


-3 


1 9 90 v 


10" 

1U 


-3 


4 


8.12 x 


io- 


-3 


—A 1 "\ * 


1U 


-3 


1 9 9^ v 




-3 


6 


8.22 x 


10" 


-3 


—A OQ y 




-3 


19 Q1 y 


1 0" 

1U 


-3 


8 


8.30 x 


10" 


-3 


—A Of, x 

L ±. \J\J A 


1 0" 


-3 


19 36 x 


10" 


-3 


10 


8.39 x 


10" 


-3 


-4.03 x 


10" 


-3 


12.42 x 


10" 


-3 


12 


8.48 x 


10" 


-3 


-4.00 x 


10" 


-3 


12.48 x 


10" 


-3 


14 


8.56 x 


10" 


-3 


-3.97 x 


10" 


-3 


12.53 x 


10" 


-3 


16 


8.65 x 


10" 


-3 


-3.94 x 


10" 


-3 


12.59 x 


10" 


-3 


18 


8.74 x 


10" 


-3 


-3.91 x 


10" 


-3 


12.65 x 


10" 


-3 


20 


8.82 x 


10" 


-3 


-3.88 x 


10" 


-3 


12.70 x 


10" 


-3 



TABLE IV: Photon inclusive corrections for p % — > irirj (see definition in Eq. (19)) as a function of 
the photon energy ujq. The masses of charged and neutral p mesons were fixed to 775 MeV's. 



